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INVESTIGATION WITH SATELLITE DATA - II 
TEMPERATURE RETRIEVALS - II 

InbrodUGtlon 

The 1976 Final Report (Fritz, 1976), hereafter designated FR-76^ and 
a published report (Fritz, 1977) described a method for retrieving atmos- 
pheric temperatures, with the aid of satellite radiance measurements and 
also nearby radiosonde, (R/S), measurements* The method was applied to 
simulated satellite radiances at ^'stations*' along a line between two 
radiosonde stations in a severe storm situation, For that purpose, the 
six CO 2 ’^radiances” for VTPR spectral frequencies were computed over each 
simulated R/S temperature sounding (VTPR - Vertical Temperature Profile 
Radiometer), The method developed in FR-76, derived a set of coeffi- 
cients, which when multiplied by the measured radiances, yielded smaller 
temperature retrieval errors than the minimum-information retrieval method. 
At the two test R/S stations, the error was in fact forced to zero. 

The purpose of this research program is to obtain accurate retrievals 
from the forthcoming VAS radiance measurements from the geostationary sat- 
ellite (VAS = VISSR Atmospheric Sounder) , Therefore similar computations 
were made with the seven VAS spectral frequencies. The results for the 
^•stations^' along the same line between the same two R/S stations used in 
FR-76 were essentially the same as for the VTPR frequencies. 

But suppose that a severe storm situation arises in an area which 
does not lie exactly along a line between two radiosonde stations; instead 
suppose the storm area is surrounded by several radiosonde stations. Also 
suppose, as is often the case, that all measurements in the six or seven 
satellite radiance channels are not really independent; instead perhaps 
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only three are independent. \'/hat would he the best method of treating such 
a situation? In that case, let's assume, for example, that four radiosonde 
stations surround the storm area and three Independent satellite radiance 
channels are available. Let one radiosonde station serve as a "first 
guess" estimate of temperature. It would then be possible to derive the 
required coefficients by solving three equations in three unknowns, These 
coefficients would achieve zero temperature retrieval error at the three 
radiosonde stations. The coefficients could then be applied to points 
insi ie or near the area surrounded by the R/S stations, As will be shown 
later, when this was done with simulated data, the errors were relatively 
small. 

It therefore seemed desirable to apply the method to real data. 

Data from VAS are not yet available. Therefore, to minimize the cloud 
contamination problem, microwave data from the NIMBUS 6 SCAMS instrument 
were used (SCA>IS = Scanning Microwave Radiometer). As might be expected, 
the errors with real data were somewhat larger than the errors with 
simulated data, but still seem smaller chan errors ordinarily obtained 
with the minimum-information method. However problems arose, related to 
the non-independence of the three SCAMS radiance channels; SCitMS often 
does not contain even three independent radiances. These problems are 
described later in this report. 

Retrievals VTPR vs VAS Spectral Frequencies 

The results reported in FR-76 were based on VTPR spectral frequencies. 
To be more applicable to the forthcoming VAS, the method proposed in FR-76 
should be applied to VAS radiance frequencies. This was done with trans- 
mittances kindly supplied by Mr. H, Fleming (NESS/NOAA) . The frequencies 
used in the VTPR and VAS retrievals are as follows: 
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Table of Frequencies (cm 


VAS 


680.0 

692.0 

703.0 

715.0 

745.0 

760.0 

775.0 

VTPR 

i 

668.5 

677.9 

695.3 

708.6 

725.5 

7iI7.6 

- 



As might have been expected, the retrieval temperature errors with VAS 
frequencies were not much different from the VTPR errors reported in FR~76. 
Tills is indicated in Table 1. If anything, the VTPR results were slightly 
better than the VAS results; e.g., for the average temperature in the layer 
1000-800 mb, the errors were 0.3K vs 0.8K, 

TABLE 1. 

Retrieval Temperature Errors : VAS vs VTPR 
Average Absolute Error (°K) 


Pressure Adjustment Method First-Guess Min-info 
(mb) Method 



VAS 

VTPR 



1000-800 

0.8 

0.3 

3.0 

2.9 

800-600 

1,0 

1.2 

7.4 

2.3 

600-400 

1.4 

1.0 

7.1 

1.6 

400-200 

1.3 

1.2 

4.0 

3.8 

200-5 

0.3 

0.2 

4.6 

0.1 


For comparison, the average absolute difference between the First Guess 
and the true temperatures, and also the average absolute error from the mini 
mum information method, are repeated from FR-76. As noted in FR-76, the er- 
rors in the proposed method are smaller than the errors in the minimum in- 
formation method, especially in the lower layers of the atmosphere. 

Retrievals N-Eguations in N-unknowns , 

As mentioned in the Introduction, a potential severe storm area will 
in general not lie along a line between two radiosonde stations, but might 











rather be surrounded by a few stations, Moreover, generally the radiances 
from the six VTPR or seven VAS COg channels are not independent. We might 
therefore have as many radiosonde stations surrounding the storm area as 
there are independent satellite radiance channels. From these we need to 
derive a set of coefficients with which to multiply the radiances in order 
to retrieve the atmospheric temperatures. In that event, we can obtain 
the coefficients exactly by solving N-equatlons in N-unknowns 5 for this, 

N equals the number of radiance channels used, 

Let us review briefly the linear retrieval method. 

The linear problem is reduced to finding a set of coefficients which 
are to be multiplied by corresponding radiances. Tlie problem is to find 
the "best" set of coefficients F(k,j) in the following relationship: 


AB (k) = EF(k,j)*AR(d) 
J 

(1) 

B^(k) = Bg(k) + AB^(k) 

(2) 

AR(j) - R^(j) - Rg(j) 

0 ) 


In these equations 

B(k) = Planck function at the pressure level, k, at a reference fre- 
quency, usually taken at 700 cm when the 15 pm CO 2 band is used in the 
satellite measurements. 

R(j) is a radiance at the frequency, j. 

The subscripts denote the following: 

"g" * a guess (or initial estimate) 
m = the measured value 
r = the retrieved value. 

Because B(k) depends only on the temperature, T(k), it is easy to 
compute T(k) from the inverse of the Planck function, once B(k) is known. 
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In 3?R-76 there wore six radiances; therefore there were also six coef- 
ficients, F(k,j) In the rainlmum-lnfortnatlon retrieval method. The method 
In FR-76 consisted of modifying two or more coefficients (at each pressure 
level, k) obtained from the minlmum-lnformation method. The new set of 
coefficients, G(k,j) yielded zero error at the radiosonde stations. \^hen 
these new coefficients were used with simulated satellite data at "stations" 
between the R/S stations, the retrieved temperatures were more accurate than 
retrievals from the minimum-information method; this was already shown In 
Table 1, above. 

If wo use only N spectral Intervals and N radiosonde stetdlons, we can 
write Eq, (1) separately for each R/S station, substituting the true value 
of the Planck function, B , for the retrieval value, B , 

C IT 

Thus for R/S station 1 we would have 

ABj.i(k) = IF(k,j)*ARj^(j) - F(k,l)’ARj^(l)H +F(k.>N) -AR^(N) 

for station 2 

= EF(k,j)-AR2(j) 

and so forth until the number of stations equals the number of frequencies 
used, N. 

Since the F(k,j) remain the same at every station, we arrive at the N 
equations in the N unknown values of F(k,j), j = 1,N, 

Test for Two R/S Stations . To test the selection of coefficients by solv- 
ing N-equatlons in N-unknowns, we used the same two stations as in FR-76; 
namely, Omaha, Nebraska and Springfield, 111. And we used two VAS spectral 
radiances; namely channels 2 and 5, or wave numbers 692 cm ^ and 745 cm ^ 
respectively. In this case, the simulated "station" at "Omaha +150 Ion" 
was used as the First Guess station, and the coefficients were derived so 


6 


as to give zero error at Omaha and at Springfield, For the four remoining 
Stations located at 50 knii 100 km, 200 km and 250 km from Omaha, the Average 
Absolute Erj'ors of the temperature retrievals are shown in Table 2. 


TABLE 2 

Average Absolute Error (°K) 



FG-True 

Retrieval 
2 channels 

3 channels 

1000-800 

0.6 

0.1 

0.2 

800-600 

2.4 

1.1 

0.5 

600-400 

2.0 

0,4 

0.3 

400-200 

1.0 

Bjjjjl 

1.9 

200-5 

1*5 


0.1 


The results obtained when only two channels were used reduced the error 

when compared with the first guess errors. Still the average error in the 

layer 800-600 mbs was more than 1°K. In an attempt to reduce that error, 

three channels were used namely, channels 2, 5, 7; channel 7 corresponds 

to wave number 775 cm”^. And in addition to the same R/S stations as were 

used before, an additional condition was imposed; namely EF(k,j) “ 1, The 

i 

justification for such a condition was discussed in FR-76; its use will, 
however, Improve the retrievals only when the coefficients are reduced in 
magnitude. In this case, it did reduce the average absolute error in the 
layer 800-600 mbs from I.IK to 0.5K. But the error in the layer 400-200 mb 
was Increased, Therefore, no conclusion about the efficacy of the added 
condition can be drawn. More exp erimence would be required before its 
merit can be evaluated. 
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Recrleval8-"~ThyGG Radiosonde Stations > Not In a Lino, 

Ab mentioned earlier, it is desirable to examine the retrievals based 
on R/S stations which do not lie along a line. The R/S stations could per- 
haps surround an area over which a storm might develop. 

We therefore selected Monett, Ho,, and Bismarck, N,D, as the test sta- 
tions, Peoria, 111, served as the "first-guess" station. Tlie temperatures 
at standord pressure levels were the average monthly temperatures tabulated 
in the climatological data for April 1973 (NOAA, 1973). The simulated sat- 
ellite "radiances" were computed for each station. With two test stations, 
we could utilize only two satellite radiance channels. In order to utilize 

a third channel, the condition BP(U,j) « 1 was also employed. The three 

3 

VAS channels used were channels 2, 5, 7 which correspond to wave numbers 
692, 7A5, 775 cm ^ respectively. 

With that combination of R/S stations and satellite radiances, we had 
three equations from which we solved for three coefficients at each level, 
Then to test the coefficients, a set of additional stations arranged in a 
circular rings were used. Again the temperatures were obtained from NOAA, 
1973. Those stations and the retrieval errors are indicated in Table 3. 

In Table 3, we note that although the First Guess error was 4,2K at 
St. Cloud, Minn., the retrieval error was only 0,7K for the layer 1000-800 mb 
Similarly, the retrieval errors were everywhere IK or less in Table 3. A 
trial with these stations for the minimum-information method, showed quite 
large errors. But even if the minimum-information error could be reduced 
to the 3K or so which is typical for that method in the lower atmospheric 
layers, the retrieval method described here would give better results, 
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Table 3 

TEMPERATURE ERROR °K 


FIRST GUESS-TRUE 

Ap 

(wb) 

PEORIA 

MONETT 

BISMARCK 

N. PUTTE 

R/iPID CITY 

DODGE CITY 

ST. CLOUD 

1000-800 

0 

-2,0 

+4,9 

2.5 

2.9 

0.2 

4.2 

800-600 

0 


3.5 

0,3 

1,3 

-2,3 

2,6 

600-400 

0 

WM 

3.6 

1,0 

2,5 

-0,7 

2,6 

400-200 

0 


2.1 

1.1 

2.2 

-0.5 

2,0 

200-5 

0 


-1.1 

1 

-0,4 

-0.5 

-0,1 

-0,8 



3 Frequencies (Channels 2, 5, 7) 


PEORIA = First Guess Station 

MONETT and BISMARCK “Adjustment Stations 
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An attempt wao qIbo made to uuc four channels; the channels used wore 
at wave-numbers 692, 701, 765, and 775. Here Peoria v;aa sRain ujed ns the 
"First-Guess" station, and, Monott, Bismarck, and North Platte, Neb,, were 
used as the test otaflons. To use four frequoncios with only these three 
test stations on additional condition was required, One possible condi- 
tion is to make f;F(k,j) •=> 1 at each level, k. This was done with the 

j 

results sho\m in Table 4(b). l/e note that, unfortunately, several errors 
were larger than I'^K, In the layer 1000-800 mb, errors of 1,7 and -1,0 
appear; in the layer 400-200 mbs one error reached 2.7°K, 

These errors were associated with large values of F(k,J), In on 
attempt to reduce the values of F(k,j), another condition was substituted 
for 2P(k,j) = 1. Instead we used 2F^(k,j) « minimum real constant, at 

1 J 

each level, k. The idea here is that if each F(k,j), is a vector com- 

2 

ponent, then the composite vector would be f(F) (k,j). By making 

4 

2 ** 

2(F) (k,J) small, we should reduce the values of F(k,j), This was in- 

j 

deed the case as shown in Table 4(d), For example in the lower levels, 

at p o 901 mb. F for v ■=> 703 was more than 6 timet? larger for the case 

?^hen 2F « 1; and for p = 299 mbs (k « 70), F was about 3 times larger, 

j 

And we further note in Table 4(c) that the errors are in fact reduced 
when the values of F(k,j) are reduced; in Table 4(c) f-h». errors are 
mainly smaller than in Table 4(b), 

Real Data 

All the results discussed above Involved simulated data. And 
those results suggest that the "adjustment method" can produce tempera- 
ture retrievals with relatively small errors, It therefore seemed de- 
sirable to proceed to real data. Unfortunately real sounding data 
from geostationary satellites are not yet available. Therefore data 
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TABLE 4. TEMPERATURE ERROR (°K) 

FOUR CHANNI5LS 

THREE STATIONS PLUS (b) EFCk,.1) “ (c) EF^Ck.j) “ min coustanL 

1 J 


(mb) 


PEORIA MONETT BISMARCK N» PLATTE RAPID CITY DODGE CITY ST, Ci %’D 


(a) FIRST GUESS-TRUE 


1000-800 

0 

1 

-2.0 

i. 

4.9 

2.5 

2.9 

0.2 

to 

800-600 

0 

-2.5 

3.5 

0.3 

1.4 

-2.3 

2.6 

600-400 

0 

-1.7 

3.6 

1.0 

2.5 i 

1 

-0.7 

2.6 

400-200 

0 

-1.4 

2.1 

1.1 

2.2 

-0.5 

2.0 

200-5 

0 

0.7 

-l.l 

-0.4 

-0,5 

-0.1 

-0.8 


(fa) ADJUSTMENT MEIHOD 


1000-800 

0 

0 

0 

0 

1.7 

0.5 

-1.0 

800-600 

0 

0 

0 

0 

0.3 

-0.5 

-0.1 

600*400 

0 

0 

0 

0 

-0.5 


0.9 

400-200 

0 

0 

0 

0 

2.7 


-1.7 

200-5 

0 i 

1 

j 

1 

0 

0 

-2.7 


1.2 

1 


(c) ADJUSIMENT METHOD 


1000-800 

0 

0 

0 

0 

1 

0.1 

■gH 

300-600 

0 

0 

0 

0 

1 

o 

-0.7 1 


600-400 

0 

0 

0 

0 

+0.6 

-0.6 1 


400-200 

0 

0 

0 

0 

-0.5 

0.5 


200-r, 

0 

0 

0 

0 

0.4 

-0.4 



(d) Values of F(k,j) 

(a) SF(k,j) = 1; (b) l’F^(k,j) « niinlraura 

i 



V “ Wave number (cm ~^) 



k 

Ap (mb) 

692 

703 

745 

775 


1 

97 

70 

901 

299 

(a) 17.1 

(b) 6.5 

(a) 2.1 

(b) 0.8 

L — 

-29.3 

4.7 

-3,3 

1.0 

1 

22.2 

4.3 

3.6 

1.3 

i 

i 

-9.1 1 

-1.8 i 

-1.4 

-0.5 
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from tho polar-orbltlnn satollitc* Nimbus 6, wore used. Furthermore In order 
to minimize the complications introduced by clouds, the microwave data from 
the SCAMS (Scanning Microwave Spectrometer) were used. 

The SCAMS data and radiosonde data for February 7, 1976 were kindly 
Supplied by Dr. N, Grody of NOAA/NESS. Data over Europe were used in order 
to minimize the time difference between Nimbus data (noon-midnight orbit) 
and the 12 z radiosonde data. The data were examined in two areas, one over 
Eastern Europe and one over Western Europe, 

Figure 1 shows an infrared clouu (and/or snow) picture, taken by sat- 
ellite NOAA 4, over Europe on February 7, 1976. The two rectangular areas, 
enclosed by solid lines, ar^i the Eastern and Western European areas for 
which retrievals were calculated in this paper. The area enclosed by the 
dashed line approximates the area of the synoptic analysis map of Figure 2. 
Figure 1 shows that some clouds were present over Western Europe; however 
near SON, lOE the sky seemed clear, -ver Eastern Europe, there also were 
some clouds, but these may have been lower than those over Western Europe, 

Figure 2 shows the synoptic situation as represented by the heights 
and temperatures of the 850 mb surface. A high contour area was located 
near SUN, 20E. This was associated with cold air further east and a 
fairly steep E-W gradient of temperature there. Over Western Europe the 
air was warmer; and more importantly, a warm area was centered near SON, 

6E so that the temperature gradient there was small.. 

Figure 3 shows the weighting functions of the three SCAMS channels. 

The two lower channels "see" the ground ; this is especially true for 
channel 3 (52,85 GHZ). We note also the substantial variation with 
nadir angle; this points up the possible error introduced by Inexact 
correction of the observed radiance for nadir angle. 
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Figures 4 and 5 show the radiance distribution from fiCAMS for Eastern 
Europe and Western Europe respectively, The patterns of temperature (Fig, 2) 
and radiance (Figs., 4 and 5) are fairly similar although the radiance pat- 
terns show more detail. This might be expected, if only because of the 
high spatial resolution of the radiance data. Over Eastern Europe, the. 
north-south and east-west gradients are evident in both radiance and 
temperature data. Over Western Europe, the gradients are more irregular, 

Tt is however interesting to note that even relatively small features ap- 
pear in both the temperature and radiance distributions. Thus the low 
radiances near 47N, lOE correspond to the tongue of cold air extending 
from Eastern Europe Into Western Europe in Fig. 2, However in the radi- 
ance data this cold area seems to bi* cut off from the main cold area; in 
Fig. 2, the spatial resolution is inadequate to show such a detail. 


TABLE 5 

Average Absolute Temperature Error (°K) 
Eastern Europe 


SCAMS (Nimbus 6) 


SIMULATED IR 

-- 

Ap 

(mb) 

F' “True 

Bi 

i 

FG-True 

Min-Info 

Adjustment 

Method 

1000-800 

5.7 

1.7 

/ 

/ 

/ 

/ 

/ 

/ 

/! 


2.9 

0.3 

800-600 

4.4 

1.5 

HI 

2.3 

1.2 

600-400 

4.4 

1.4 


1.6 

1.0 

400-200 

2.8 

1.3 

4.0 

3.8 

1.2 

















There are many radiosonde stations ovcir Europe, And since the SCAMS 
makes measurements at three frequencies designed for temperature retrievals 
(Staelln, et al. , 1975), three radiosonde stations in the middle of the 
SCAMS data over Europe were selected. This was done separately for Eastern 
Europe and separately for Western Europe. In each region one additional 
station was selected as "first-guess" station. We therefore have three 
equations like Eq. (1). [For the microwave region it is appropriate to 
use AT instead of AB, since, to a close approximation, B is a linear 
function of T, at a given wavelength.] 

From the three equations, three coefficients F(k,j) were computed 
at each level, k. To test the method, the SC.\MS radiance measurements 
over a set of additional radiosonde stations were used. With these 
radiances, and using the coefficients, F(k,j) just derived, the tem- 
peratures were computed (or retrieved) at the additional radiosonde 
stations. The retrieved temperatures, T^, were then compared with the 

corresponding radiosonde temperatures, to determine the errors in T . 

r 

Eastern Europe The results are shown in Table 5 for Eastern Europe 
(Fritz, 1977(b)), Table 5 shows the average absolute temperature error. 

In the layer 1000-800 mb, where large errors are often made in the mini- 
mum informatior method, the error was 1.7K; the average absolute dif- 
ference between the "First-Guess" temperatures and the true temperature 
was 5.7K. For comparison. Table 5 shows also the errors from earlier 
reports (Fritz, 1976, 1977). We note here again that even with real 
data the errors seem smaller than are usually encountered with the 
minimum information method. 

The errors for Individual stations are shown in Fig. 7, along- 
side the line for channel 3. The errors were at most 2°K except for 
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the two stations which were close to the satellite horizon. For those 
the error was 2.1 and 2.6°K. The reason for the somewhat larger error 
In the stations near the horizon may be In the corrections applied to 
the radiance. According to Dr. Grody (verbal communication) a correc- 
tion is made to the radiances for nadir angle of the area viewed from 
the satellite. This correction is an a priori correction, based on a 
statistical sample from around the world obtained prior to the satel- 
lite launch. Obviously this correction can only be approximate and 
would be more inaccurate for large nadir angles than for small nadir 
angles . 

However, the errors with real data were larger than the errors 
encountered with s'*jnuiated data when the adjustment method was used. 
This is, of course, to be expected; with real data new sources of 
error appear. Not only are there errors in the radiance data, (which 
were, to some extent, also introduced in simulated data) but there are 
errors in the "truth" the radiosonde temperatures contain errors. 

Some error is introduced by inevitable time and space mismatch 
of the satellite radiance and radiosonde data. Clouds may influence 
the radiances somewhat, although in the microwave data, the cloud 
effect is minimized because the clouds are mainly transparent to 
microwave radiation. The ground emissivity may be different from 
unity and may vary from station to station. Therefore, with real 
data, "errors" [which are really differences from the estimated 
"true" values] should be larger than for simulated data. 

Still the retrieved temperature errors in Table 5 are small 
even for the real data. 
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It should ho pointed out that even for Eastern Europe, some trial and 
error was required in order to reduce the errors; that is, the selection 
of the first guess station as well as the test stotions from which the 
coefficients F(k,j) are derived must be carefully selected. 

This can best be explained if we recall that the solution for the 
three values of F(k,j) Involve a fraction; thus for F(k,j) corresponding 
to channel 3, 


F(k,3) 


AR^l 

AR,2 

AT^ 

AR21 

AR22 

AT2 

AR31 

"^32 

AT3 

“u 


AR ^_3 

^^21 

AR22 

AR23 

AR3, 

"^32 

AR33 


Num 

Den 


Eq, (6) 


In order to arrive at small values of F(k,j), — a necessity if temperature 
errors are to be small — it is desirable to have the numerator small and the 
denominator large. If the AR's for each channel are highly correlated, the 
denominator will approach zero, and F(k,j) will tend to be large. It is 
therefore desirable to select a set of stations for which the AR's are 
independent and at the same time the AT' s are relatively small. 

Western Europe 

Unfortunately, for Western Europe, the method of using three equations 
in three unknowns did not work well. This was caused by the fact that the 
three SCAMS radiances were not sufficiently independent over Western Europe. 


l-Jben the radiances are not independent, the coefficients, become 

large, as was shown above. With large F(k,j) the errors become large. 

Tlie lack of independence of the SCAMS data over Western Europe is 
shown in Fig, 6. Figyre 6 is a plot of AR, for channel 3, against AR for 
channel 2, and for channel 1, The figure shows that the radiances are too 
highly correlated (correlation coefficient a, 0,9). This leads to large val- 
ues of F(k,j) for all sets of test stations tried. As a consequence, we 
have essentially only one independent channel and the adjustment method 
requiring 2 or more channels will not work. 

Regression Method . 

With only one channel it may still be possible to obtain retrievals 
with small errors by regression methods. This would be possible if the 
temperatures were correlated with the radiances. And, indeed, the values 
of ATj^Q_g[= T^Q_g - wevQ well correlated with the corresponding 

values of AR for channel 3 and channel 2, Such a relationship is illus- 
trated in Fig. 7 for Eastern Europe, The correlation between AT and AR 
are high for both channels 2 and 3. For channel 1, the correlation was 
poorer. 

Similar curves are available for Western Europe, There the correla- 
tion is not as high largely because the variability of the temperatures 
was less than in Eastern Europe. 

Because of the good correlation between AT and AR, we should expect 
reasonably accurate retrievals, even from one radiance channel. Therefore 
a regression method was applied to the radiances to obtain the temperatures . 

This was done with a statistical package available in the University 
of Maryland computers, namely, Subroutine RLSEP, which selects a regression 
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model uaing a forward stopwloe olgorithra, computes the coefficients In the 
rcgrespUm equation and olso various other statistics. In RLSEP it is pos' 
sible to force all the radiances into the multiple regression; RLSEP can 
also be used so that only statistically "significant" radiances are used 
in the regression. Both methods were employed. 

Western Europe 

It turned out that only one radiance channel was significant for 
Western S>irope. Usually, the results with one channel alone, were about 
as good as the result when all three channels were used. This is lndi~ 
cated in Table 6. From Table 6, we note that the percent variance ex- 
plained in the dependent data is 84% when only one radiance channel was 
used. This was improved to only 88% when all three channels were used, 
More importantly when the regression equations shown in Table 6 were ap- 
plied to independent data, the Average Absolute Error (AAE) was about 
1.8K whether one radiance channel or three channels were used. In fact 
with the three channels, the error, if anything, was even slightly larger. 

Table 6 also show that when all stations, the dependent together with 
those which had been used as independent stations, were used to obtain re- 
gression equations, the conclusion was the same; namely, the addition of 
two channels to channel 2 added almost nothing to the percent variance 
explained or to the standard deviation of the residuals. 

Eastern Europe 

Even though the method involving three equations was successful 
over Eastern Europe, the regression was employed there also. Table 6 
shows the regression equations. The percent variance explained was very 
high, 97%; and the standard deviation of the residuals was small; namely, 
about 1.3K. The regression subroutine selected both channels 1 and 3 


TABLE 6, Regression Equations, ^T^gOO-SOO Radiance; ?; Var = Percent Variance Explained 

S.D. = Standard Deviation of Residuals, AAE = Average Absolute Error 
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as significant for Eastern Europe, Apparently channel 2 did not add signi- 
ficantly to the explanation of variance. This might have been expected 
also from Fig. 7, which suggests that AR(3) and AU(2) are highly correlated, 
because both ore highly correlated with AT, The fact that the regression 
equations over Western Europe are different from those over Eastern Europe 
was also observed by Grody and Pellegrino (1977). 

Discussion . 

The regression method when used with only one radiance channel depends 
on the good relationship between AT and AR. iflien one considers that T for 
the layer 1000-800 mbs contributes very little to the radiance in channel 3, 
the question arises as to why highly correlated with AR(3) , 

the A radiance in channel 3, and with AR(2). Such a correlation probably 
requires that the air temperature over a deep layer should also be corre- 
lated with 

To test this the average temperature for the layer 1000-800 mb, T-- o> 

XO— o 

was plotted against the average temperature in the layer 800-600 mbs. The 
result is shown in Fig. 8, for Eastern Europe. As expected, the relation- 
ship is quite good. Thus, when ^^Q_g Is cold, 5g__g is cold also; this of 
course generally agrees with meteorological experience also. Major cold 
air outbreaks involve cold air over much of the troposphere. Also lati- 
tudinal temperature variation in winters has at least the same sign over 
most of the troposphere. 

As a consequence, not only is AR(3) related to is 

also highly correlated with AR(2), 

Regression vs Adjustment Method. The regression method when applied over 
a limited geographical area can give as good results as any other method. 
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But tho regression method requires a relatively large number of radiosonde 
Stations in order to make the regression coefficients stable. Moreover 
using relationships from some other orea or some other time leads to in- 
creased errors (Grody and Pellegrino, 1977). Therefore, for the times and 
place for which retrievals are desired it may be necessary to acquire R/S 
data over a relatively large area. This in turn may moke the synoptic 
situation too variable to achieve good results over a small severe storm 
area. Therefore, where possible it may still be desirable to use the ad- 
justment method, with three or four equations involving three of four 
radiosondes around the area under study. ^Vlicn this is not possible, the 
regression method can be used but at the expense of Involving data from 
a larger area. 
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Fig. 1. Infrared picture from satellite NOAA-4. Rectangles 
with solid lines are areas over Eastern Europe and 
V^cstern Europe for which temperature retrievals were 
made with SCAMS date. Irregular dashed line bounds area 
of synoptic chart in Fig. 2. (Courtesy Grody and 
Pellegrino) . 
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2. Synoptic 850 mb chart over Europe. Solid lines are contours 
dashed lines are isotherms. Feb. 7, 1976» 12Z. 
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Fig. 3. Weighting functions for SCAMS channels in the 
oxygen band. The three channels were designated 1, 

2, 3; channel 1 corresponding to the uppermost channel, 
35.45 GHZ, and channel 3 to the lowermost, channel 
52.85 GHZ. (Courtesy Grody and Pellegrino) 








minus 200. Number *’50” corresponds to 230K. Dashed lines 
bound the area %rith relatively small satellite nadir angle 
The black dots are the locations of radiosonde stations. 
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SCAMS NIMBUS- 6 FEB., 7, 1976 (~I0= 20 GMT) 
WESTERN EUROPE (~42-58N,l-l5E ) 



Relation of radiances in channel 3 with radiances in 
channel 2 and in channel 1. 








EASTERN EUROPE 



oetween average temperature in layer 1000-800 mb with 
average temperature in layer 800-600 mbs. Line drawn by eye 






